Coherent spin control by electrical manipulation of the magnetic anisotropy 
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High-spin paramagnetic manganese defects in polar piezoelectric zinc oxide exhibit a simple almost 
axial anisotropy and phase coherence times of the order of a millisecond at low temperatures. The 
anisotropy energy is tunable using an externally applied electric field. This can be used to control 
electrically the phase of spin superpositions and to drive spin transitions with resonant microwave 
electric fields. 
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Couplings between magnetic and electric degrees of 
freedom give rise to such fundamental phenomena in con- 
densed matter as multiferroelectricity [TJ [2], unconven- 
tional superconductivity [3] and spin-density- waves [4], 
and are key to proposed future technologies such as high 
sensitivity metrology [5 and quantum-dot-based quan- 
tum information processing (QIP) [6J. The control of spin 
states using electric fields rather than magnetic fields is 
particularly valuable for QIP [7], because electric fields 
can be applied on short length scales, and down to the 
scale of qubit separations in a device [8HTD] . In this Letter 
we identify a class of electrically-controllable spin qubits: 
high-spin magnetic defects in polar semiconductor hosts. 
In one member of this class, the manganese substitutional 
defect in a crystalline zinc oxide host, we demonstrate co- 
herent oscillations of the spin state driven by d.c. electric 
field pulses and spin transitions driven by resonant mi- 
crowave electric fields. 

Electrical control has been considered in various can- 
didate physical systems for spin qubits, both experimen- 
tally and theoretically. In lithographically defined GaAs 
quantum dots, spin-orbit coupling permits coherent con- 
trol of an electron spin by applying a resonant high fre- 
quency voltage to one of the confining gates [11]. Stark 
shifts in spin- splittings, originating from hybridization 
with excited state orbit als, have been observed in self- 
assembled quantum dots [12 and in bound donors in 
semiconductors [13] [14] and in defects in diamond [15] , 
while frustrated spin triangles exhibit spin-electric cou- 
plings via modulation of exchange interactions [16] [17] . 
Here, our approach is to manipulate the electrostatic 
environment of a high-spin paramagnetic defect, man- 
ganese, in a polar piezoelectric host material, zinc oxide, 
thereby controlling electrically the spin Hamiltonian of 
the defect. 

The spin Hamiltonian of manganese defects in zinc 
oxide (ZnOiMn) is well-known from continuous- wave 
ESR [18 . The polar environment lifts the spin degen- 
eracy of the S = 5/2 Mn 2+ ions which are substitutional 
on Zn 2+ sites, while the 7 = 5/2 nuclear spin of 55 Mn is 
coupled to the electron spin by the hyperfine interaction. 
The Hamiltonian, including a magnetic field which lifts 



the remaining Kramer's degeneracies is 

H = fji B S • g • B + g N n N i B + S A I + H ZFS (1) 

where the axial anisotropy, or zero field splitting, term 

Hzfs = -DSl (2) 

Higher order effects, such as fourth order electron spin 
anisotropy and nuclear quadrupole terms, have been de- 
tected [18] , but they are small and not relevant for the 
following discussion. The resulting energy spectrum is 
shown as a function of magnetic field applied at a small 
angle to the anisotropy axis in Fig. [TJa). The absorp- 
tion in response to an oscillating magnetic field of fre- 
quency 9.7 GHz, applied perpendicular to the static field, 
is shown in Fig.[TJb); the peaks are magnetic dipole elec- 
tron spin resonance transitions corresponding to the se- 
lection rule Am s = ±1. 

The spin states are highly coherent. Fig. 1(c) shows 
the relaxation times for electron spin population (Xi e ) 
and coherence (T2 e ) as a function of temperature, mea- 
sured on the m s = ±1/2 transition (□ symbols in Fig lb). 
The relaxation times increase rapidly below room tem- 
perature with Ti e ~ T~ 3 and T 2 e ~ T -2 , saturating at 
~ 0.2 s and T 2e « 0.8 ms below 10 K. The cubic tem- 
perature dependence of T\ e is consistent with scattering 
by a population of phonons, an indication that lattice 
distortions couple to the Mn spin. 

The total Mn density is 4.0 ± 0.1 x 10 16 cm" 3 , corre- 
sponding to a mean separation of 17 nm between near- 
est neighbours. The thermal spin populations lead to 
an effective dilution of on-resonance spins by a factor of 
(2S + 1)(2J + 1) = 36, so spin centres in the same state 
are separated on average by 58 nm. The quadratic T 2e 
temperature dependence is characteristic of decoherence 
driven by spin flips (i.e. T ie relaxation) of off-resonant 
neighbouring spins [19] [20], and the saturation below 
10 K is consistent with magnetic dipole couplings be- 
tween more distant spins in the same sub-level. 

The relaxation times are independent of the nuclear 
spin projection to within experimental errors, but as ex- 
pected for a phonon-driven mechanism, T 2e depends on 
the m s states involved in the superposition. The coher- 
ence times for m s = ±1/2 : ±3/2 superpositions and 
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FIG. 1. (Color online) (a) The energy level structure of the system as a function of magnetic field, (b) The ESR spectrum 
comprises five groups (Am s = ±1) of six (21+1) lines. Weakly allowed (Ara s ,Araj) = (±1,±1) transitions are visible 
near the centre of the spectrum, c) Electron spin lifetimes T± e and T2 e for the m s — —1/2 : 1/2 transition. Fits are to 
Tie = AT 2 /(exp(-E/k B T) + 1) modelling spin relaxation via phonons, T2 e ~ (Ti e ) 2//3 modelling relaxation caused by spin 
flips of neighbouring Mn spins. 



m s — ±3/2 : ±5/2 are, respectively, about 0.65 and 0.58 
times that for m s = —1/2 : ±1/2. The ratios are not 
strongly temperature-dependent. 

The large moment of the 6 S5/2 electronic configuration 
allows us to perform Rabi oscillations with periods down 
to ~ 5 ns (depending on the m s states involved) with 
standard apparatus, allowing us to perform ~ 4 x 10 6 
single qubit operations within the relaxation time of this 
system. 

The origin of the axial anisotropy of the Mn spins, 
parameterised by D in Eqn. [I] is the electric field experi- 
enced by the Mn defects in the polar ZnO. Applying an 
additional external electric field E along the polar (100) 
axis allows us to modify the zero-field-splitting parameter 
such that D = Dq ± k,E, and it is possible to manipu- 
late the spin via this externally-controllable parameter. 
Fig|2|a) shows spin state control using a d.c. electric field 
pulse embedded into a Hahn echo sequence [13] [21] . The 
initial tt/2 pulse generates a superposition between ad- 
jacent m s levels; the electric field adjusts their splitting 
for a time r p , causing a phase 

Acj) = olkEt v (3) 

to accumulate in the superposition; the 7r-pulse refocuses 
static inhomogeneities; and an echo forms whose phase 
reflects the manipulation applied by the electric field 
pulse. Phase shifts in quantum spin coherences, gener- 
ated by sweeping the d.c. electric field E between and 
4 x 10 5 V/m for a pulse duration of r p = 6 /is, are con- 
verted to populations and are shown in Fig. 2(b) for the 
five coherences between adjacent m s states. The electric- 
field-induced shifts in energy are proportional to m 2 s1 so 



the m s = ±5/2 : ±3/2 superpositions respond most 
strongly to E (a = 4 in Eqn. [3|, the m s = ±3/2 : ±1/2 
respond more weakly (a = 2), and the ±1/2 superpo- 
sition is unaffected (a = 0) to first order. We find 
only a first order phase shift from the applied electric 
field [13]. The rate of accumulation of phase gives a di- 
rect measurement of the coupling to the electric field, 
k = 20.8 ± 0.2 rads-VVm" 1 . 

The oscillations in Fig. [2^b) constitute coherent ma- 
nipulations of the Mn spin using an electric field. The 
magnitude of the effect that we observe suggests that the 
externally applied electric field does not simply add to the 
internal electric field of the ZnO. Assuming that the zero- 
field splitting in the absence of an applied electric field, 
Do = 1.585±0.006x 10 9 rads -1 , originates from an inter- 
nal electric field of a typical magnitude, ~ 10 10 Vm _1 , 
the magnitude of k that we measure demonstrates that 
the external applied electric field is several orders of mag- 
nitude more effective at modifying D than the internal 
electric field. This is probably because the dielectric and 
piezoelectric response of the ZnO contributes to the mod- 
ification of the local environment of the Mn ion, thereby 
enhancing the effect of the applied electric field on the 
spin anisotropy. While demonstrated here for Mn de- 
fects in ZnO, this effect should be observable for any 
S > 1 spin in a lattice whose crystal-field is modified by 
an external electric field. 

The capability to tune spins into and out of reso- 
nance with a cavity using an electric field is highly desir- 
able for information processing applications [8]. Elec- 
tric fields of order 10 6 Vm _1 , that might be applied 
by microscopic lithographically-defined structures, would 
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FIG. 2. (Color online) (a) A Hahn echo sequence enclosing 
a voltage pulse to generate phase shifts of the spin coherence, 
(b) The phase of the recorded spin echo shifts when an electric 
field impulse is applied. For the m s < transitions the elec- 
tric field generates a positive phase shift (i.e. the quadrature 
channel Q, plotted in green, lags the in-phase channel, /, red), 
whereas for m s > a negative phase shift occurs, (c) The 
recorded oscillations are unwrapped using 
and plotted against the electric field impulse. The phase shift 
is proportional to the electric field impulse, indicating that 
the shift in D is linear in E. 

generate shifts in the spin resonance frequencies of or- 
der 10 MHz. Comparing this tuning bandwidth to the 
~ 100 kHz linewidths of contemporary superconducting 
microwave resonators [22 , we can envisage hybrid sys- 
tems in which d.c. electrical tuning via an array of sur- 
face gates selectively brings chosen spin qubits into and 
out of resonance with a superconducting cavity bus. 

We have seen so far how d.c. electric fields can shift 
energy levels and cause phases to accumulate in preex- 
isting spin state superpositions. However, it is also possi- 
ble to create superpositions and excite coherences using 
resonant a.c. electric fields. For a magnetic field applied 
along z at an angle 9 to the anisotropy axis, the zero field 
splitting Hamiltonian becomes 

#zfs(0) = -^costf + ^sintf) 2 . (4) 

If D is modulated at high frequency by applying an os- 
cillatory electric field along the anisotropy axis, i^zFS 
can generate resonant transitions between spin states m s . 
Fig. [3^ a) shows the transition strengths driven by an 
electric field of frequency 9.7 GHz as a function of mag- 
netic field strength and orientation 0. When the mag- 
netic field is perpendicular to the anisotropy axis and 
the oscillating electric field (0 = 90°), H ZFS = -DSl = 
— ■^D(S+ + S-) 2 . Hzfs contains pairs of raising and low- 
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FIG. 3. (Color online) (a) (Theoretical) Electric dipole spin 
transition probability per unit time (indicated by colour den- 
sity) in ZnO:Mn as a function of the angle between the 
static magnetic field and the anisotropy axis for an oscillatory 
electric field of frequency 9.7 GHz. Single quantum coher- 
ences are plotted in black, double quantum coherences in red. 
(b) (Theoretical) The corresponding conventional magnetic 
dipole electron spin resonance transition probability, for an 
oscillating magnetic field of the same frequency applied per- 
pendicular to the static magnetic field, (c) (Experimental) 
The measured electric dipole transition strengths for small 
angles 0. The upper inset shows the experimental geometry; 
the oscillatory electric field strength in the rectangular TE102 
mode is indicated by colour density (red is strong, blue is 
weak). The lower inset shows the raw data for the derivative 
absorption spectra as a function of magnetic field for small 
angles 0. 

ering operators, so double quantum coherences may be 
generated by a resonant oscillating electric field, i.e. tran- 
sitions can be excited for which Am s = ±2; these transi- 
tions are shown in red in Fig.[3^a). For a general non-zero 
#zfs contains first and second powers of the raising 
and lowering operators, driving both single and double 
quantum coherences (transitions for which Am s = ±1 
or ±2). For = 0, when both fields are parallel to the 
anisotropy axis, the total spin Hamiltonian is secular and 
no transitions are excited. For reference, the conventional 
electron spin resonance magnetic dipole transitions are 
shown in Fig. |3jb), for a static magnetic field applied at 
an angle to the anisotropy axis and an oscillatory mag- 
netic field applied perpendicular to the static magnetic 
field. 

Fig. j3^c) confirms that we can excite electric dipole 
transitions between spin states. The sample is posi- 
tioned on the conducting wall in an electric field antin- 
ode (amplitude ~ 500 V/cm) of a rectangular TE102 res- 
onator (Q ~ 5000), such that the oscillating electric field 
(which is perpendicular to the conducting wall) is paral- 
lel to the anisotropy axis. The orientation of the static 
magnetic field, which for = is parallel to the os- 
cillatory electric field, can be varied in our apparatus 
by up to about 5 degrees. The inset shows the deriva- 
tive absorption spectrum obtained as a function of mag- 
netic field strength as is varied; the magnetic field 
range is chosen to examine the single quantum coher- 
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ence transitions m s = +5/2 : +3/2 because, along with 
m s = —5/2 : —3/2, these should exhibit the strongest 
electric dipole transitions. For = no transitions are 
visible, and as the magnetic field is rotated away from 
the anisotropy axis, the transition strengths increase in 
proportion to the component of the oscillating electric 
field perpendicular to the magnetic field. This is exactly 
the behaviour expected for the electric dipole transitions 
as predicted in Fig. |3|a). 

We can rule out the possibility that the transitions 
in Fig. [3jc) are magnetic dipole transitions. Firstly, the 
sample is small compared to the cavity dimensions, so the 
oscillating magnetic field that it experiences is vanish- 
ingly small. Secondly, any such oscillating magnetic field 
that the sample does experience is parallel to the con- 
ducting wall and is therefore perpendicular to the static 
magnetic field for = 0; in this geometry the magnetic 
dipole transitions are maximised and, to first order, in- 
sensitive to small changes in 6. 

We have shown that the electron spins of substitutional 
manganese in zinc oxide exhibit very long coherence 
times and that their states may be manipulated rapidly 
using d.c. and a.c. electric fields. Owing to the possibility 
of applying electric fields over much shorter length scales 
than magnetic fields, such electric field control of qubits 
offers dramatic advantages for quantum information pro- 
cessor architectures [23ti25] . The spin-electric coupling is 
several orders of magnitude larger than is found in other 
electric-field-sensitive spin qubit candidates p~3j [15] , and 
large enough that electrical control of the spin state is 
possible on time scales comparable with magnetic con- 
trol: the coupling of the Mn 2+ spin to the electric field of 
a plane wave in bulk ZnO is about one order of magnitude 
smaller than the strength of its coupling to the magnetic 
field; the electric field coupling could be enhanced by 
confining the electromagnetic wave in a structure with 
a higher characteristic impedence than the bulk. The 
origin of the electric field coupling that we identify here 
should be applicable to many high-spin magnetic defects 
in polar materials. In particular, we would expect that 
electrically charged Fe 3+ impurities in ZnO should show 
even larger couplings than Mn 2+ impurities, which are 
isoelectronic with the Zn 2+ . Other polar materials, such 
as the ferroelectric strontium titanate, should also be ex- 
pected to host magnetic defects exhibiting strong electric 
field couplings. 

We thank W. Hayes, A. Heinrich, T. Lancaster and 
J. Moeller for helpful discussions. We thank J.J.L. Mor- 
ton for discussions and materials. This research was sup- 
ported by the UK EPSRC. A. A. thanks the Royal Society 
for support. 



[4] 
[5] 
[61 

[7] 

[8] 

[9] 

[TT 
[12 

[13 
[14 

[15 

m 
w 

[18 
[19 

[20 



[21 
[22 

[23 
[24 
[25 



(2007) 

G. Gruner, Density waves in solids (Addison Wesley, 
1994) 



J. Zhai, Z. Xing, S. Dong, J. Li, and D. Viehland, Applied 
Physics L etters 88, 062510 (2006) 

J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 

A. Yacoby, M. D. Lukin, C . M. Marcus, M. P. Hanson, 

and A. C. Gossard, [Science 309, 2180 (2005) 

L. C. Bassett, F. J. Heremans, C. G. Yale, B. B. Buckley, 

and D. D. Awschalom, |Phys. Rev. Lett.| 107, 266403 (Dec 

2011) 

A. Andre, D. DeMille, J. M. Doyle, M. D. Luki n, S. E 
Maxwell, P. Rabl, J. Schoelkopf, and P. Zoller, |Nature| 
Physics 2, 636 (2006) 

H. Bernien, L. Childress, L. Robledo, M. Markham, 
D. Twitchen, and R. Hanson, Phys. R ev. Lett . 108, 
043604 (Jan 2012) 

S. T. Ochsenbein and D. R. Gamelin, Nature Nanotech 



nology 6, 112 (2011) 
K. C. Nowack, F. H. L. Koppens, Y. V. Nazarov, and 
L. M. K. Vandersypen, Science) 318, 1430 (2007) 
J. M. Smith, P. A. Dalgarno, R. J. Warburton, A. O. 
Govorov, K. Karrai, B. D. Gerardot, and P. M. Petroff, 
|Phys. Rev. Lett~l 94, 197402 (May 2005) 
F. R. Bradbury, A. M. Tyryshkin, G. Sabouret, J. Bokor, 
T. Schenkel, and S. A. Lyon, Phys. Rev. Lett7| 97, 176404 
(Oct 2006) 

A. Morello, J. J. Pla, F. A. Zwanenburg, K. W. Chan, 
K. Y. Tan, H. Huebl, M. Mottonen, C. D. Nugroho, 
C. Yang, J. A. van Donkelaar, A. D. C. Alves, D. N. 
Jamieson, C. C. Escott, L. C. L. Hollenberg, R. G. Clark, 
and A. S. Dzurak, Nature 467, 687 (2010) 



E. V. Oort and M. Glasbeek, Chemical Physics Letters 
168, 529 (1990), ISSN 0009-2614 



M. Trif, F. Troiani, D. Stepanenko, and D. Loss, Phys. 
Rev. Lett.| 101, 217201 (Nov 2008) 



M. T rif, F. Troiani, D. Stepanenko, and D. Loss, Phys. 
Rev. B|82, 045429 (Jul 2010) 



[1] D. Khomskii, Physics 2, 20 (Mar 2009) 

[2] N. A. Spaldin and M. Fiebig, [Scie^ 309, 391 (2005) 

[3] P Monthoux, P. D., and L. G.G., [Nature] 450, 1177 



A. Hausmann and H. Huppertz, J. Phys. Chem. Solids 
29, 1369 (1968) 

A. Barbon, M. Brustolon, A. Lisa Maniero, M. Ro- 
manelli, and L.-C. Brunei, Phys. Chem. Chem. PhysTj l, 
4015 (1999) 

A. M. Tyryshkin, S. Tojo, J. J. L. Morton, H. Riemann, 
N. V. Abrosimov, P. Becker, H.-J. Pohl, T. Schenkel, 
M. L. W. Thewalt, K. M. Itoh, and S. A. Lyon, Nature 
Materials! 11, 143 (2012) 

W. Mims, The Linear Electric Field Effect in Paramag- 
netic Resonance (Clarendon Press, Oxford, 1976) 
P. J. Leek, M. Baur, J. M. Fink, R. Bianchett i, L. StefTen, 
S. Filipp, and A. Wallraff, |Phys. Rev. Lett! 104, 100504 
(Mar 2010) 

H. Wu, R. E. George, J. H. Wesenberg, K. M0lmer, 
D. I. Schuster, R. J. Schoelkopf, K. M. Ito h, A. Arda- 
van, J. J. L. Morton, and G. A. D. Briggs, Phys. Rev. 
Tett]l05, 140503 (Sep 2010) 

Y. Kubo, F. R. Ong, P. Bertet, D. Vion, V. Jacques, 
D. Zheng, A. Dreau, J.-F. Roch, A. Auffeves, F. Jelezko, 
J. Wrachtrup, M. F. Barthe, P. Bergonzo, and D. Esteve, 
|Phys. Rev. Lett~] l05, 140502 (Sep 2010) 
X. Zhu, S. Saito, A. Kemp, K. Kakuyanagi, S. Kari- 
moto, H. Nakano, W. J. Munro, Y. Tokura, M. S. Everitt, 
K. Nemoto, M. Kasu, N. Mizuochi, and K. Semba, Na- 
ture 478, 221 (2011) 



